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ABSTRACT 



Two puzzles associated with open clusters have attracted a lot of attention - their formation, with densities and velocity 
dispersions that are not too different from those of the star forming regions in the Galaxy, given that the observed Star 
Formation Efficiencies (SFE) are low and, the mass segregation observed / inferred in some of them, at ages significantly 
less than the dynamical relaxation times in them. Gas dynamical friction has been considered before as a mechanism for 
contracting embedded stellar clusters, by dissipating their energy. This would locally raise the SFE which might then 
allow bound clusters to form. Noticing that dynamical friction is inherently capable of producing mass segregation, 
since here, the dissipation rate is proportional to the mass of the body experiencing the force, we explore further, some 
of the details and implications of such a scenario, vis-a-vis observations. Making analytical approximations, we obtain 
a boundary value for the density of a star forming clump of given mass, such that, stellar clusters born in clumps which 
have densities higher than this, could emerge bound after gas loss. For a clump of given mass and density, we find a 
critical mass such that, sub-condensations with larger masses than this could suffer significant segregation within the 
clump. 
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CN i. Introduction 

>■ . . . 

^^ Stellar clusters can be viewed out to far distances and hence are important probes for our study of the universe. Most 

l/~i fetars seem to form in clusters, and hence, understanding clusters becomes an essential requirement for understanding the 

lO formation and evolution of galaxies. Two puzzles associated with open clusters that have attracted a lot of attention are 

^""^ ■- their formation, with densities and velocity dispersions that are not too different fr om those of the star forming regions 

■^ in th e Galaxy, given that the observed Star Form ation Efficiencies (SF E) are low d Tutukov 1978; Hills 1980: Mathieu 



O [19831 ; lElmegreenI Il98^ iLada. Margulis fc Dearborn 1984 : Adams 200dt iGever fc Burkert 2001; Baumgardt fc Kroupa 
C^ 2007; Goodwin 2009) a nd, the mass segregation observed / inferred in some of them (Ide Griis et al. 20021: ISirianni et aL 
2002t iGouhermis et al.1 120041 : IChen. de Griis fc Zhaal[2007HConverse fc Stahleill2010l:lHasan fc Hasadl201lh. at ages sig^ ^ 



^ nifica ntly less than their dyna mical relaxation times (iBonnell fc Daviea 119981: iMcMillan. Vesperini fc Portegies Zwarlj 
• !-( I2007t lAlhson et aLll2009l l2010f . iMaschberger fc Clarkeir2011l: lOlczak. Spurzem fc Henning||201lD . In this paper, by mak- 
/\^ ing analytical approximations, we explore the efficacy of gas dynamical friction, operating during the embedded phase 
H of stellar clusters, in solving these two puzzles. 
. . -Star clusters form, hidden from sight, embedded in the dense, dark and cold cores of giant molecular clouds (jLada fc Ladal 
[2003). A star forming gas c loud converts less than ten percent of its mass into stars and then, in less than ten million 
years (jPalla fc StahlerJl2000D loses the gas which had bee n gravitationally binding the st ellar cluster. A cluster of stars 
will be bound if the total energy of the stars is negative (jChandrasekhar fc Elberdll972[ ). For mass loss that is instan- 
taneous, analytical approximations (|Hiiy[l98(i) show that the Star Formation Efficiency (SFE), i.e. the fraction of the 
total mass of the gas cloud that is converted into stars, should be greater than 0.5, if a bound cluster is to emerge after 
gas loss. On the other ha nd, if the mass loss is gradual, bound clusters similar to those observed could emerge for SFE's 
greater than ~ 0.3 (JMathieul [1981 rElmegreenlll983t iLada. Margulis fc Dearbornlll98l) . Note however that this value 
is still much higher than the observed values for the global SFE. The SFE in dense cores which are forming clusters 
seem to be ^ 0.2 — 0.3, while the observed values for the SFE in star forming clouds are an order of magnitude smaller 
iLadaJf Lada"2003'). 

Observations show also that, emb edded clusters can start showing mass segregation in less than a few million years 
(jSchmeia. Kumar fc Ferreirall2008r ) and that the Orion Nebula Cluster is not old e nough for the observed segregation of 
its massive stars to be due to dynamical relaxation by two body encounters alone (jBonnell fc Davieal 19981) . Simulations 
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indicate that the Pleiades cluster started out with signific ant mass seKregation in th e embedded phase itself, with high 
mass stars preferentially concentrated towards the center ([Converse fc Stahleij|2010r ). The dynamical relaxation time in 
Pleiades is twice the presen t age of the cluste r. Also, at the same time, some cluster s like Taurus and Trumpler 16 do 
not show mass segregation ([Parker et alll201 iL IWolk et alJIMllllHasan fc Hasaiill201lf ). 

Numerical investigations of the dependence of the survivability of a cluster on the Star Formation Efficiency has shown 
that, a significant fraction of the stars can remain bound, after rapid gas loss, for an SFE > 0.3 and, for slow gas 
removal, up to 50 % of the cluster can remain bound for SFE values down to 0.15 — 0.2 ("Geve r fc BurkertI 120011 : 
iBaumgardt fc Krou pa 2007; Goodwin 2009). However, if most open clust ers are remnants of e mbedded clusters that have 
- con sequent to gas removal - lost one half or more of their natal members ([Bastian fc Goodwin 2006; G oodwin fc BastianI 
[200a), the similarity in the mass functions of em bedded and open clu sters and, the similarity in the IMF's of embedded 
and open clusters become difficult to understand ([Lada fc Ladal[2003[ ) . In this context, scenarios in which bound clusters 
can form, without losing too much of their natal stellar content and, with (possibly) mass segregation become interesting. 
Gas dynamical friction due to an embedding gaseous medium has been invoke d, in understa nding many phenomena; at 
various scales - from that of planetary systems to that of clusters of galaxies (iKim fc Kim 2009). In t he cas e of stars, 



while considering mechanisms that change the momenta of stars in stellar systems, IChumak fc Chumalg ([1976[ ) conclude 
that, dynamical friction in a dust cloud would be the dominant momentum changing mechanism for a star, when there 
are no massive scattering c entres and the speed of the star does not exceed a l i mit. Making use of the Chandrasekhar 
dynamical friction formula ([ChandrasekhaHri943l ial). [Saivadpour. Deiss fc Kegel! ([1997 1 had explored gas dynamical fric- 
tion as a means to produce contraction of embedded clusters, thus raising the SFE locally, which could then lead to the 
form ai on of bound clu s ters. They had suggested also that, the same might be responsible for mass segregation ([Deiss et all 
|1998( ). iGorti fc BhattI ([19951 ) had explored the effect of gas dynamical friction on pre-stellar clumps by numerical simula- 
tions and also, it has been noticed that gas dynamical friction can reduce the time scale for dynamical mass segregation 
([Er et al.ll2009ll . However, no general analysis of such a scenario, or its details and implications vis-a-vis observations, 
has been made. Here, using analytical approximations, we explore a scenario in which gas dynamical friction could be 
playing a significant role, in the early evolution of embedded stellar clusters, as the chief mechanism which can cause 
them to contract and form more strongly bound configurations and, also cause mass segregation in them. We check the 
plausibility of the model by examining its consistency with observations. 

In section 2, using simple virial considerations, we check that, a scenario in which embedded clusters undergo con- 
traction, meets the dynamical constraints set by observations. In section 3 we consider dynamical friction in star forming 
gas clouds. By requiring that the gas dynamical friction time scale within an embedded cluster be less than the time that 
is available before gas is expelled from the cluster, we obtain a boundary value for the density of a star forming clump of 
given mass, such that, stellar clusters born in clumps which have densities higher than this, could emerge bound after gas 
loss, and for a clump of given mass and density, we find a critical mass such that, subcondensations with larger masses 
than this could suffer significant segregation within the clump. In section 4 we compare our results with observations. 
In section 5 we discuss some of the possible implications of our scenario and, their accord with observations. In section 
6 we give a brief summary of our results. 

2. The SFE as a constraint on cluster formation scenarios -analytical approximations 

We now proceed as follows, for checking, whether the scenario we propose, wherein a nascent cluster contracts within 
the parent gas cloud, meets the constraints set by observations. For a gas cloud that converts a fraction e of its mass 
into stars and then loses the remaining gas, iHillsl ([1980|) had, by simple virial considerations, obtained expressions for 
the ratios of, the radius (rg), velocity dispersion (ao) and density (po) of the gas cloud, with the radius (rf), vel ocity 
dispersion (ct/) and density (pj) of the stellar cluster, assumed virialized after gas loss (see also ([Mathieulll983f) '). as 
c/ZcTo = \fi-l \/tj jrij and Pf / po = e/(rf/rQ)^, where rf /rQ equals e/(2e — 1) for instantaneous gas loss and equals 1/e 
for adiabatic gas loss. Making a similar analysis for an embedded stellar cluster, with e < 0.5, that contracts to a radius 
r' — pro, where p the contraction factor is < 1 we get, rf/r^ = {p/{l — ((1 — e)/e)p^)) for instantaneous gas loss and 
r//ro = p(l + ((1 — e)/e)p^), if the gas is lost adiabatically. 

For the first two cases, we see from the corresponding expressions for rf/ro that, in the case of instantaneous mass 
loss, bound clusters may form only for SFE's > 0.5 and in the case of slow mass loss, bound clusters may form, techni- 
cally, for all values of e > 0. In the case where the embedded cluster contracts to within a radius r' — pro, we see that, 
even if gas loss is instantaneous, bound clusters may form, for all values of e < 0.5 also, if p is < ((e/(l — e)) . For 
P < ((^/(l ^ e)) ' , the contracted cluster will have an apparent SFE > 0.5 ([Saivadpour. Deiss fc Kegellll997t ). 

In a real situation, a cluster will be situated in a galaxy, and hence, will be able to retain only the stars that are 
within the tidal radius. Our analysis, thus may be applied to those clusters, that do not fill their tidal radii and also, 
are young enough for their dynamics to have not been affected by either internal processes, like mass loss from the stars 
or external processes like interactions with molecular clouds. The typical value for the density of a star forming region 
is ~ 10^ cm^"^ (Elmegreen 1983; Bcrgin fc Tafalla 2007). And, taking for typical cl usters, a de nsity enhancement of 50 
with respect to the density of field stars - which is of the order of 4 x 10~^^gcm~^ ([Coxlll999r ) - we get a typical value 
for Pf/po ^ 0.005. The typical value for the observed velocity dispersion for clusters is ~ 0.5kms~^ and for molecular 
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clouds it is of the order of a few kilometer per second, giving a typical value for (Jf/oQ ^ 1/4 — 1/6. 

We now estimate a typical value for the contraction factor p, by obtaining an estimate for ry/rg using observations 
as, r-f/ro — 1.7(e/0.025)^/'^((p//po)/0-005)^^/'^ and then, solving for p using the relation rf/ro = (p/(l — ((1 — e)/e)p^)) 
(in the case of instantaneous gas loss) or, rf/ro = p(l + ((1 — e)/e)p"^)) (in the case of adiabatic gas loss). The estimate 
for p may be used to calculate Ca = 1/(1 + ((1 — e)/e)p'^), the apparent SFE of the cluster, when it has contracted 
to such an extent that, if it loses the gas now, it can remain bound, with a radius rj (after virialization), such that 
rf/ro — 1-7. For rf/ro = 1-7 and e = 0.025, we get the required values for the contraction factor p as, p — 0.278 
in the case of fast gas loss and p = 0.426 in the case of slow gas loss. These values for p, correspond respectively to 
apparent SFE's ~ 0.544 and 0.25 for the embedded stellar cluster, in its most contracted phase, just prior to losing the gas. 

In our scenario, we would expect at least some embedded stellar clusters, to show such high (apparent) SFE's, when 
near to the end of their embedded phase. Stellar densities that imply SFE's upto 0.47 have bee n reported, in the dense 
embedded clusters observed in the p Ophiucus, IC348 and Orion- Trapezium star forming regions (jWilking fc Ladal (|l9 8tl) 
and references therein). We interpret these, as apparent SFE's, arising from the contraction of the particular nascent 
cluster. Thus these high SFE's that have been observed, lend support to our scenario. 

As e is varied from 0.008 to 0.08 the estimates for the appropriate apparent SFE's, as given by our model, go from 
0.54 to 0.576 for fast gas loss and, from ~ 0.2 to 0.3 for slow gas loss. Assuming that rf/ro has a larger value, will make 
the apparent SFE values that correspond to it smaller, and vice versa. As an aside we notice that, the value for the ratio 
af/ao, calculated from the observed value for the ratio of the densities as, (Jf /(Jo = l/8(e/0.025)^/'^((p//po)/0-005)"'^/^ 
is consistent with its value determined independently from direct observations - a reflection of the fact that molecular 
clouds, as well as bound clusters, are both roughly in virial equilibrium, a condition which was assumed in deriving these 
relations. 

Thus, we see that, our results linking r//ro,CT//tTo and pf/po, obtained under the assumption that the stellar cluster 
contracts within the parent gas cloud, are broadly consistent with observations. This conclusion of ours, is not very sen- 
sitive to the precise value of e used in the calculations, as long as it is low. With the above conclusion as our motivation, 
we now explore, by making analytical approximations, the efficacy of gas dynamical friction, in bringing about such a 
contraction of embedded stellar clusters. 

3. The gas dynamical friction time-scale 

Molecular clouds which are observed to be the sites of star and cluster formation, have a hierarchical structure, of dense 
cores within clumps, the clumps themselves bei ng within clouds - the density getting lower as the scale becomes larger, 
from cores to clumps to clouds (see for example (jBergin fc Tafallall2007t )). For a clump of mass ra, and velocity v, moving 
in a medium with density Pgas, the gas dynamical friction time scale tgdf is ~ {E/{dE/dt)), where dE/dt is the rate of 
energy dissipation due to gas dynamical friction and E is the kinetic energy of the clump. 

The retarding force, when a massive perturber that is interacting gravitationally moves through a fluid, may be de- 
termined by considering the i nteraction of th e body with its own gravitationally induced over-density wake. For density 
wakes with small amplitudes. lOstriken (|1999D using time-dependent linear perturbation theory, determined the drag force 
on a point mass perturber, moving in a straight line, through a uniform, infinite, gaseous medium. For a mass m, moving 
with a speed f, in a medium with density pgaa and sound speed c^, she obtained an expression for the force as 

F.,„ = -l^^^H^^v/(M) (1) 

where /(M) = l/21n((l + M)/(l - M)) - M for M < 1 and = l/2/n(l - (l/M^)) + ln{vt/rrmn) for M > 1. Here M ^ v/cs 
is the Mach number, t is the time for which the perturber has been moving in the medium, and rmin is a minimum 
radius introduced to avoid singularity in the force evaluation. Later workers have relaxed the constraint on the perturber 
being point-like, to a perturber which has no boundary but is merging into the surrounding medium; that of a linear 
trajectory to considering circular orbits; and that of a uniform medium to one that is radially stratified (Kim & Kim 
2009 and references therein). These authors found that equation (1) is generally applicable with appropriate changes to 
the factor [vt/rmin)- 

Thus we get 

—- = fcFLi„.v = k "-—^ ^/ M , 2 

dt V 

where k is a numerical factor that accounts for departures from the linear theory. For a clump, moving within a 
cloud of total mass Mc and radius R, with a speed v, parametrized by the virial speed for the cloud as v — fivmriai, 
putting pgas — SgasP whcrc, the gas fraction egas — Mgas/Mc, Mgas being the mass in the interclump medium and p is 
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Afc/(47r/3)i?3, we get 

'"f ekG^megasMJiUy ^ > 

Real cores, clumps and clouds are not spherically symmetric, have density gradients, and also su bstructure at a variety 
of sca les. However, potentials are rounder and smoother than the underlying density distribution ([Baumgardt &: Kroupal 
I2007D . and relaxing the assumption of sphericity, or homogeneity, introduces only changes by factors of order unity in th e 
various terms of the virial equation (iChandrase khar fc Elbertlll972t ISom Sunder fc Kochhaijll985l : I VerschuerenI 1 1 990( ) . 
Hence, for the term Vyi^iai , we use the expression for the virial speed for a smooth homogeneous sphere of mass M^. and 
radius R viz. ^{2>/b){GMc/R), and write 

''"^ - 20V 57rfce,,,/(M)'^™^^ " ^-"^^fce^^JlM) '^''°^^' ^^> 

where tcross — {Gp)^^'"^ is the crossing time in the cloud. 

We rewrite the expression for the gas dynamical friction time-scale (Eqn.4) as 

'^gdf — '^f^^cross y'^ ) 

where 

0.038/33 



kegasfiM) 
and the mass ratio factor /i, is given by 



(6) 






m 



3.1. Production of bound clusters and the occurrence of mass segregation 

The chief effect of gas dynamical friction is to decelerate objects to sonic speeds. In molecular clouds, turbulence has 
been observed to be dominated b y large modes, w hereas sound speeds are rather low and subvirial (see for example 
iKirk. Johnstone fc Tafallal (J2007[ ); lKirk et aj (J2010h '). We may expect that a cluster will be bound after gas loss, if typi- 
cal stars - i.e. those around the peak of the I MF, which c arry most of the mass in a cluster - can lose energy in a time 
shorter than the gas expulsion time-scale tE (|GieleJl2010 f). In our case this would ensure that all higher mass objects 
would also be bound. 

Thus, in our scenario, with tcross '^ {C'p)~^ , we get the condition for the formation of a bound cluster as 

15Tfl{pMQpc-3)~^ <tE{Myr)- (8) 

This yields a critical density for the parent gas cloud, such that the cluster emerges bound after gas loss, as 

P > Perineal = 225(T^/t_E(M2;r))^M0pC~3. (9) 

For given p, r and Mc, we may rewrite the above inequality to give a critical mass rricritieai such that, all objects with 

m > mcritieal = 15(TAfc/t£(Myr) ) (PMopc-3 )" " (10) 

would be significantly slowed down by gas dynamical friction. 

3.1.1. Estimates for the various factors 

For a pre-cluster, we take i^, the length of time f or the pre-stellar embedded phase (iLeisawitz. Bash fc Thaddeuslll989f ). 



when gas dynamical friction would be operative (jSanchez-Salcedo fc Brandenburd 2001), as 5Myr. Th is is less than a 



quarter of the pre- main sequence contraction time for a one solar mass star ( Charbonnel et al]ll999[ ). Stellar winds, 
which can push away the ambient gas, thus reducing gas dynamical friction, are initiated when the stars reach the main 
sequence. Mole cular outflow s which have the same effect, have dynamical time scales that are much smaller, of the order 
of O.lMyr only (|Franklll999l ). 

We may expect a bound cluster to form if cores with mass es down to the orde r of 1 M(n , which c an produce stars 
with masses u p to the peak o f the IMF (Larsen 1982; Weidner fc Kroup a 2006; Wei dner. Kroupa fc Bo nnel 201Qi), which 
is « 1/3 Mq (jChabrieiJ I2003D . have tgdf < 5Myr. This is because they will now have ample time to get decelerated to 
sub- virial speeds by gas dynamical friction, even as they migrate to the minimum of the gravitational potential. 
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We now make an estimate for r, for star forming cores in clumps, and for clumps in molecular clouds. We first con- 
sider /(M). The observed density and velocity structures in mole cular clouds are consistent with supersonic turbulence 
driven at the scale of the cloud itself (jOssenkopf fc McLowll2002[ ). In determining the Mach number M for a perturber 
moving through a turbulent medium, the speed of small scale turbulence may b e used instead of the sound speed , as the 
speed wi t h wh ich pressure perturbations are transmited through the medium (jSaivadpour. Deiss &: Kegellll997[ ). From 
lOstrikeii (119991 ) we see that, from ^ M'^/S for M << 1, /(M) rises sharply to ~ ln{vt/rmin) — 2 across M = 1 and tends 
to ln{vt/rmin) for M >> 1. We see that t he time scale will be shorter for motion that is supersonic, since /(M) is > 1 
only for supersonic motion ()Ostrikeilll999D . It will be shorter for objects that have comparatively smaller speeds, since a 
smaller value for (3, implies, both a smaller amount of energy that need be dissipated, and a smaller separation between 
the perturber and its wake. It depends on the density of the parent cloud via the crossing time, and is shorter for denser 
clouds. It decreases as the gas fraction in the cloud increases. Also, it depends directly on the mass of the parent cloud - 
the amount of energy that has to be dissipated will be proportional to Mc, since the virial speed is proportional to vMc 
- and inversely on the mass of the perturber, which decides the dissipation rate. It depends inversely on /(M), which is 
related to the speed with which pressure forces can redistribute the density enhancement in the wake. 

We consider now a distribution of speeds for the particles. Those objects which have supersonic speeds with respect to 
the large scale streaming motion of the gas surrounding it, can be significantly slowed down, by gas dynamical friction. 
For subsonic objects, the gas dynamical friction forc e is depressed due to the fact that pressure forces are able to restore 
the density distribution about the perturber, sooner (I Ostrikenll999f ). We see that the chief effect of gas dynamical friction 
would be to slow supersonic objects down to sonic speeds. Super-virial objects are likely to escape from the system. So, 
from here onwards, we put /? = 1 in the expression for tgdf- 

We now consider fc, the correction factor for non-linearity. iKim fc KimI (|2009f ) investigating the gas dynamical fric- 
tion force in the non-linear regime, by high resolution hydrodynamic simulations, had measured the strength of the 
induced gravitational perturbations, by a dimensionless parameter A = (rB/rg) = (Gm/c^rs) (tb is the Bondi radius). 
A, which corresponds roughly to the ratio of the perturbed density at a distance r^ from the perturber, to the back- 
ground density, is >> 1 in the non- linear regime. Following them, in our case, we determine the correction necessitated 
by non-linearity, as follows. For a typical protostellar object of mass IM0, taking r^ ~ lOOAU -the radius of a protostellar 
nebula, and with Cg ^ lkms~^, which is the order of the speed of sound/ random motions in molecular clouds, we get 

A = 900(kms-Vcs)2(m/Mo)(AU/rs) (11) 

to be ~ 9. iKim fc KimI ()2009f ) obtain fc in terms of a non-linearity parameter ry = {A/{yP — 1)) as, fc — (F/Fim) = 
(r//2)~*''*^ where F is the gas dynamical friction force in the non-linear regime. We convert the limiting values they 
obtain for ry, into limits on the Mach number M for the protostellar object in a star forming gas cloud, as follows. For 
M < 1 and for M > 3.7 (i.e. r] < 0.7) {F/Fun) - 1 and for 3.7 > M > 2.35 (i.e. 0.7 < r; < 2), (F/Fun) -- 1-6 ~ 1.0. 
For 2.35 > M > 1.04 (i.e. 2 < 77 < 100), {F/FLm) ^ I ~ 1/5. We see that for cores which are supersonic, but not 
greatly so, non-linearity introduces a reduction. For this range of speeds, we get a typical value for this reduction factor 
by evaluating it at M = 2 as fc = (77/2)""-^^ ^ 0.83. For clumps, which are not compact like cores - with typical masses 
~ IOOMq and sizes ~ few light-year - we are in the linear regime, and fc = 1. Thus we have fc ~ 0.83 and €gas ~ 0.9 
for cores in clumps and, fc = 1 and e^as ^ 0.75 for clumps embedded in molecular clouds. Here the gas fraction within 
clouds has been taken as the complement of the mass fraction for clumps in clouds (|Lada fc Ladall2003t iKauflmann et all 
[2009). 

For a perturber, which is moving in a turbulent medium, the wake would not be fully developed. Taking ln{vt/rmin) ~ 3, 
i.e. for a size of the wake ~ 20 times the size of the perturbe r , /(M ) ^ 1, and we find that r is '-^ 0.05. This may 
be compared with the value l/iir ~ 0.08 obtained bv lOstrikerl (|l999t) . as the value for the coefficient relating the gas 
dynamical friction time scale to, the mass ratio and the typical time scale in the system, while considering the decay of 
the near-circular orbit of a massive perturber, moving in a massive, spheri cal cloud with a singular isothermal sph ere 
profile, which result, was itself found consistent with the results obtained bv lSanchez-Salcedo fc Brandenburg (J200lt ). in 
three dimensional simulations, of a gravitational perturber moving in a gaseous sphere. 

4. Comparison with observations 

The formula we have obtained above, for the critical density PcriUcai, such that, if the cloud gas density p is > PcriUcah 
gas dynamical friction would be significant, for clumps of mass m, moving in a cloud of mass M^ — m/ p,, may be used 
to investigate the degree to which gas dynamical friction is significant for cores in clumps and for the clumps themselves, 
embedded in molecular clouds. 

For a typical clump of mass 500 M© - which, given a SFE ~ 0.1 (jLada fc Ladal l2003l ) . could be the progenitor of a 
typical cluster of mass 50 M^:^ - we get a critical d ensity ~ 1.2 x lO^cm"^. This is an order of magnitude higher than the 
density of typical clumps (jBergin fc Tafallall2007l) , which have densities ^ 10'^ — 10"* cm~'^ . Also we get the critical mass in 
typical clouds which have masses 10^ — 10''' M0 and densities 50 — 500 cm^'' as, 30 < nicr < 10'^ M0 , the lower value being 
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for less massive clouds with a higher dens ity and vice vers a. Cores which have m > rricr, may lose significant amounts 
of energy, and could congregate (see also iKirk et all ()2010 f)') at the potential minima, raising the apparent SFE there. 
The compact cluster that forms, can then emerge bound after gas expulsion, if the local SFE, gets raised sufHciently 
(see discussion before). Thus clusters which were born in sufficiently dense clumps, could emerge bound, and also show 
mass segregation (since protostars having higher masses and hence, comparatively lower values for /x, would be more 
strongly affected by gas dynamical friction), while those born in clumps with comparatively lower density, would form 
associations. For a clump of mass 10^ Mq, Pcriticai is ~ 5 x 10^ cm~'^; i.e. an order of magnitude larger, than densities in 
the upper range of values for typical clumps. 

As a particular example, we now compare the predictions of our theory, with observations of the nearest and well studied 
Orion Nebula Cluster, which is young enough for the mass segregation observed in it to h ave been very fast JR eggiani ct aj 
20111). It is ob served that only stars with masses > 5M0 are segregated in the ONC (jHillenbrand fc Hartmannlll998l) . 



Huff fc Stable r (2006) estimate the mass of the parent cloud of the ONC as 6700 M©, within a radius 2.5 pc, and we 



get the critical mass (Eqn 9), for the parent cloud of the ONC as 100 M©. From the correlation observed in embedded 
clusters, between the total mass of the stars in the cluster and the mass of the most massive star in it, the most m assive 
star a clump of mass 100 M0 can produce is » SMfn (Wcidner fc Kroupall2006l : IWeidner. Kroupa fc Bonne]|f2010f ). if we 
assume a star formation efficiency of 0.3 ( Lada fc Ladal 12003V This value is interestingly close to that for the lowest 
mass observed to be segregated in the ONC. It may be noted that the above mentioned correlation, seen between the 
total mass of the stars in a cluster and the mass of the most massive star in it, need not conflict with the observations, 
either of isolated O stars or, with numerical simula tions whic h suggest that, clusters form by random sampling from a 
universal IMF with a fixed stellar upper-mass limit (|Lamb et a l. 2010). In the first inst ance, these cou l d be ' run aways', 
and in the second instance, the correlation would still hold, albeit statistically (see also Oh fc Kroupal (J2012t )). 
The density of the parent cloud of the ONC is only ~ 2 x lO'^cm^'^. which in our scenario, is less than the critical density 
(< 10^ cm"-^) required for a bound cluster to form, from a cloud of mass 6700 M0 (Eqn 10 ). Thus, in our scena rio, the 
ONC will evolve into an unbound OB Association and this is consistent with the predictions iHuff fc Stabler ('2006') make. 



regar ding the fate of the ONC, from observations of the velocity dispersion of the stars in the cluster ([Jones fc Walken 
MM)- 

From the expression for the critical mass (Eqn 10), we see that, in our scenario, the lowest mass which will be seg- 
regated, depends directly on the total mass of the clump, inversely on the square root of the density of the clump and, 
will be smaller, for smaller speeds of random motions of the gas. This last dependence, on the level of small scale turbu- 
lence, comes from the fact that gas dynamical friction can slow condensations down to the sound speed (here the speed of 
small scale turbulence), and the lower the speed of the condensations the smaller will be the t wo body relaxation time (see 
Eqn 12 below). The above expectations from our theory, are bourne out by observations (jSchmeia. Kumar fc Ferreiral 
l2008t IChavarria et al.ll2010l ). wherein, the most centrally condensed and mass segregated clusters are found in clouds 
with the lowest Mach n umber and vice versa. By making a comparison between the age of the cluster and the migration 
time scale for the stars. IChavarria et al.l (|2010l ) attribute, the mass segregation present in the embedded clusters observed 
by them, to gas dynamical friction. The densities of these star forming clumps (which are > lO^cm"'^), are less than the 
critical densities (Eqn 9) estimated for their masses (which are > 5 x 10"^ M©), and we find that, within them, stars with 
masses ^ 1 Mp, d o not show any concentration towards the center. We do not expect any of the f our clumps obse r ved b y 
IChavarria et al.l ()2010[ ) to produce bound clusters. However the clump DR21(0H) observed bv .Bontemps et aLl (|2011f ). 
with a mass of 7000 Mq and a radius of 0.3 pc has a critical mass of only 4.2 Mq and, we consider it highly probable that 
this cloud will eventually produce a bound cluster. 

5. Discussion 

In our scenario, clusters with large masses, could form at the centers of sufficiently dense molecular clouds, by the 
congregation of the more massive clumps in it, due to decelaration by gas dynamical friction. Observationally some 
molecular clouds as well as Giant Molecular Clou ds, do show, such dense, cluster-bearing concentrations of gas in them 
(jWilking fc Ladal[l983HMaiz-Apellanidl200l[ ). and lWilking fc Ladal (|l983t l comment on their pro bable significance in the 
formation of bound stellar clusters. Also, gas that harbors clusters, is seen to be more den se ("Kauff mann et aLll2010t) 
and more highly clumped or substructured than gas that does not bear a cluster (|Lada fc L ada 200iJ). In the present 
scenario, the cluster will bear the imprint of the merged clumps and cores, as substructure in its stellar distribution, till 
it completely relaxes. This conforms wit h what is observatio nally seen, since very young, (< IMyr) star clusters, seem 
to show significant levels of substructure (jLada fc Ladal[20p3f) and, observations show that cores and young clusters have 
substructure and are cool (sub-virial) (see for example (jSchmeia. Kumar fc Ferreiral 120081 : iKirk. Johnstone fc Tafallal 
[2007)). This is consistent also, with the observed difference, between the fractal dimensions of young clusters with inter- 
nal substructure (that has still not been erased by dynamical relaxation) and, that of molecular gas, and we note that, a 
more clustered distribution of the densest star forming gas, within the parent molecular clouds - as e xpected in our sce- 
nario - has already been suggested as a probable reason for the above difference in fractal dimensions (jAlfaro fc Sanchea 
I2OIII) . Mergers, need not erase any mass segregation that was originally present within the clumps, since dynamical 
simulations have shown that, in the merging of clust ers, which are themselves mass segregated, th e mass segregation 
will be retained, rather than get erased in the merger ([McMillan. Vesperini fc Portegies Zwarlll2008l ). In this context, it 
may be interesting to note that the comparatively less dense Taurus complex, does not show any central concentration 
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or mass segregation. 

The critical densities above which we exp ect the formation of b ound clusters / mass segregation, are consistent with the 
inter-clump medium densities, for which iGorti fc BhattI (|l995f ) notice mass segregation in their numeri cal simulations 
studying the effect of gas dynamical friction on prestellar clumps and, ISaivadpour. Deiss fc Kegel! (|1997l) obtain bound 
clusters in their analysis. The larger densities, which in our scenario are a prerequisite for the formation of bound 
clusters, are also consistent with the general conclusion that the production of a bound cluster must require special 
physical conditions, to acco unt for their rarity; only about 4 — 7% of all embedded clusters seem to give rise to bound 
clusters (|Lada fc Ladall2003( ). Also lKauffmann et al.l ()2010f ) report that, between clouds which have the same radii, those 
which harbor clusters are more massive than those which do not show a concentration of stars. They also report that 
while clouds do seem to follow a general mass -size relation of the form Mc{r) = 4OOM0rpc"'^'^, there are outliers whose 
masses (and hence densities) are an order of m agnitude larger. In our scenario, such denser clouds could give rise to 
bound clusters. Goddard. Bastian fc KennicuttI (|2010; ) in a study, on the fraction of clusters that escape infant mortality 
find that, in the galaxies they studied, this fraction is strongly proportional to the Star Formation Rate (SFR) density 
of the parent galaxy. Since the SFR density is related to the surface gas density by the Kennicutt-Schmidt law, they 
conclude that the survival fraction of clusters is strongly dependent on the surface gas density. Since we expect clouds, 
and hence clumps, in regions of larger surface gas densities to be denser, this is consistent with the expectation from our 
scenario. Our model is consistent also with the ob served co r relati on, that the two categories of massive clusters - bound 
or 'starburst' and unbound or 'leaky' - noticed bv iPfalzned (|2009[ ) obey; that bound clusters are fou nd in regions o f high 
density viz. the Galactic Center or in the spiral arms, and unbound ones in the lower density areas (|PfalzneiJl2011|) . She 
even suggests, as a hint for the origin of the two categories, this obvious difference in their locations. It is also interesting 
to note that while, among these clusters, the young, bound ones have stellar densities > lO^cm"^, the initial stellar 
densities for the young, unbou nd ones are < 5 x lO'^cm"'^. The high initial densities quoted for the bound 'starburst' 
clusters - ~ 10^ — 10^ cm~'^ - ()PfalzneiJl2009t iPfalzner fc Eckartll2009[ ) are consistent with the densities expected for the 
compact clusters, that in our scenario, should form in clouds with densities greater than the critical densities we obtain 
(Eqn 9) (see also Maiz-ApeUaniz (2001)). 

The gas dynamical friction time scale is inversely proportional to the mass of the body which is being decelerated 
by it. Hence it is possible that, relative segregation, between the stars / cores / clumps of various masses may occur 
and, some signature of this might survive, even after relaxation, in the clusters evolving out of dense gas clouds. We 
note that, the br own dwarfs associated with a cluster are distribu ted in a region mu ch larger than what is occupied by 
the cluster itself (jLuhmanI 120061 : iKumar fc Schmeial 12007; Kouwe nhoven^ Brown fc Kaper 2007), and clusters do seem 
to have associated with them a coro na of low mass s tars ffSharm a et all 120061 ; lLuhmanii2006i : iKumar fc Schmeial l2007f . 



iKouwenhoven. Brown fc KaDeill2007); see also Gorti fc Bhatt (199g)). Such a 'core-halo' structure has been noticed for 



some Massive Young Clusters too (jMaiz-Apellanizj 120011) 



It is seen from numerical simulations that, clusters that ar e dynamically cool at birth can remain bound, down to 
SFE's - 0.05, if they lose the gas while they are stiU cool (|Goodwinl [19971 l2009t ISmith et aOl201lh and that , clus- 



ters that are cool and/or have substructure, can mass segrega te in a very short time sc ale (lAllison et al.l 20091 2010: 
iMaschberger fc Clarkd l2011t iOlczak. Spurzem fc Hennind 120111 ) . However, the study bv iPelupessv fc Portegies Zward 
(J2012r ). which extended these studies by including the potential due to the gas, found that, only a remnant with an IMF 
that was possibly skewed, rather than a mass segregated cluster, is left after gas expulsion. That bound clusters may 
be obtained, by starting with 'cool' initial conditions, and losing the gas before the stars acquire too high speeds due 
to collapse, may be expected from the fact that, between being born with a very low velocity dispersion, and acquiring 
virial speeds by collapsing within the gravitational potential of the gas cloud - within about a dynamical time, which is 
^ half a million year, for a density ^ lO'^cm"'^ - (jLada. Margulis fc Dearbornlll984l )). the stars can, for a short time, 
satisfy the condition for binding - that the total energy of the stars, taken alone, be negative. However, this would imply 
that the survivability of clusters would be proportional to the dynamical time, which is inversely proportional to the 
square root of the density. This would produce an anti-correlation between the cluster survivability and the gas density, 
which is counter to observations, as explained above (Goddard. Bastian fc Kennicutt 2010). 

The scenarios mentioned above, which produce fast mass segregation, may be understood by noticing the following. 
The formula for the two b ody dynamical relaxation time for a stellar cluster, whose stars are moving with virial speeds 
(jBinnev fc Tremainell2008l ). has to be modified when we consider a cluster that is embedded in gas. With e equal to the 
ratio of the total mass in stars to the total mass of stars plus gas we get 

^Erelax — Tvi T7 Across v^^j 

SmaA' 

where a = /3^/e. e is < 1 for an embedded cluster and, is equal to 1 for a cluster free of gas. The relaxation time would be 
smaller for - a) larger SFE /gas free clusters b) sub- virial speeds for the stars i.e. /3 < 1, and also c) within sub-structure, 
since in this case, the number of stars as well as the crossing time would be much smaller compared to those for the 
cluster as a whole. These conclusions are in concurrence with the results of the simulations mentioned above, wherein, 
the signature of dynamical relaxation viz. mass segregation, is reported to occur fast, in simulations which start with 
suitable (as per our analysis above) initial conditions. Our model, explores the possibility of gas dynamical friction being 
the probable cause of the occurrence of 'suitable' initial conditions, and also delineates the situations in which such initial 
conditions might arise. 
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In our scenario, the difficulty posed by the large dynamical relaxation times in embedded clusters, is avoided by the fact 
that, gas dynamical friction is capable of poducing mass segregation explicitly - due to its mass dependence, as well 
as implicitly too - by reducing the dynamical relaxation time, through its ability to decelerate condensations to sonic 
speeds, which in star forming clouds, are subvirial (i.e. by making /? go < 1) (Eqn 12). The efficacy of gas dynamical 
friction may be enhanced, under the following conditions. In turbulent molecular clouds, the decay time for small scale 
turbulence, will be smaller than that for large scale motions, if nit < 3, where Et(kt) oc /ct"™* is the energy density of 
turbulent motions per unit wave number range, kt being the wave number. For example, for a Kolmogorov spectrum 
of turbulence, mt — 5/3. Hence the sound speed, which is here taken as the typical speed of small scale turbulence in 
the cloud, can decrease over time and go over to the actual speed of sound for quiescent gas, in situations where the 
turbulence is neither forced nor fed, and clumps and cores, which are slowed down to sonic speeds by gas dynamical 
friction, can be quite sub-virial. This will reduce the relaxation time. The dying down of turbulence can cause the 
friction to increase also, by allowing longer wakes to develop. Also, for a cloud with a power law density gradient ex r"™"* 
(m^ < 3), between a small inner radius within which the density is a constant, and a large outer radius beyond which 
the density may be considered negligible, we get Mc{r) < p{r) >~^"(x r^~^~ , < p{r) > being the mean density within 
a radius r. Thus the gas dynamical friction time scale is likely to be smaller than the global mean value towards the 
centres of real clumps and clouds, as the product McP~^''^, which occurs within the equation for tgdf (Eqn. 4), decreases 
with r. 



Though gas dynamical friction can caus e decay of binary orbits (iGorti fc Bhafrg 1996a) and, binaries and stellar clusters 



can have s ignificant effects on each other (|KrouDall995tlKrouDa. Petr fc McCaughreanll999l : lDe lgado-D onate. Clarke fc Batd 
l2003l : [Parker. Goodwin fc Allison 2011r). we do not consider the possible role, that either binaries or tidal fields (see for 
example iBallesteros-Paredes et all ()2009[) ). could play in the formation of bound clusters, as they l ie beyond the sc ope 
of this paper. The model avoids the need, to take recourse to density dependent IMF's or SFE's (jKruiissenll2012f) . in 
accounting for various aspects of stellar clusters. A comparison with such models will be made in a future work. 
The model offers the following interesting possibility for the formation of Globular Clusters - by the contraction and 
coalescence, in a similar scenario, of dense cores, within clumps, within clouds, within very massive, very dense mo lecular 
clouds, that perhaps formed via turbulent fragmentation during galaxy formation (see also iMaiz- Apellania (J200lh ). 

6. Summary 

Gas dynamical friction can decelerate masses and also explicitly as well as implicitly, promote mass segregation. Com- 
paring the gas dynamical friction time-scale for sub-structures within star forming clouds, with the time available before 
gas expulsion takes place from an embedded stellar cluster, we obtain a boundary value for the density of a star forming 
clump of given mass, such that, stellar clusters born in clumps which have densities higher than this, could emerge 
bound after gas loss, due to the locally elevated SEE. For a clump of given mass and density, we find a critical mass such 
that, subcondensations with larger masses than this could suffer significant segregation within the clump. We compare 
our results with observations and also other work in the field, and also discuss the implications of our scenario for the 
formation of bound clusters, vis-a-vis observations. 
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